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We investigate the magnetic response in the quantized spin Hall (SH) phase of a layered- 
honeycomb lattice system with intrinsic spin-orbit coupling Aso and on-site Hubbard U . The 
response is characterized by the parameter g = 4Ua 2 d/3, where a and d are the lattice constant 
and interlayer distance, respectively. When g < (cr^/i) -1 , where a% y is the quantized spin Hall 
conductivity and fj, is the magnetic permeability, the magnetic field inside the sample oscillates 
spatially. The oscillation vanishes in the non-interacting limit U — >• 0. When g > (a^fi)" 1 , the 
system shows perfect diamagnetism, i.e., the Meissner effect occurs. We find that a superlattice 
structure with large a is favorable for observing these phenomena. We also point out that, 
as a result of Zeeman coupling, the topologically protected helical edge states show weak 
diamagnetism that is independent of g. 

KEYWORDS: magnetic response, quantized spin Hall effect, Kane-Mele model, intrinsic spin-orbit interac- 
tion, electron correlation, layered-honeycomb structure, superlattice, topological BF term, 
superconductivity 



1. Introduction 

The aim of this paper is to discuss the magnetic re- 
sponse of the quantized spin Hall (SH) phase of corre- 
lated electrons on the layered-honeycomb lattice. 

It has been shown that the quantum SH effect occurs 
in a non-interacting electron system with intrinsic spin- 
orbit coupling Aso: like the Kane-Mele (KM) model. 1-3 ' 
The investigation of the quantum SH system with on-site 
Hubbard U is now one of the current topics, 4-6 -* and the 
phase diagram has been obtained. 5 ' 6 ' In ref. 7, magnetic 
response in the quantized SH phase of the layered KM 
model with U (spin-conserving limit of the topological 
band insulating phase in ref. 6) was discussed. Such a 
model can be applied to correlated electrons in the sys- 
tem with honeycomb layers, such as some Ir-based ox- 
ides. 4 ' The correlation is characterized by the parameter 
g oc Ua 2 d in the low-energy long-wavelength regime (o; 
lattice constant, d; interlayer distance), and the London 
equation for the Meissner effect has been obtained in the 
large g limit. 7 ' The discussion on general g has not yet 
been given. The role of the topologically protected he- 
lical edge states, which is the hallmark of the quantum 
SH effect, 1,2,8-10 ' has also not been taken into account. 

In this study, we clarify the magnetic response for gen- 
eral g in the quantized SH phase. First, we take into ac- 
count contributions from bulk states and later consider 
those from the helical edge states. When g < (cr^ 2 ^) -1 , 
where a s xy and \i are the quantized spin Hall conductivity 
(SHC) and magnetic permeability, respectively, the mag- 
netic field inside the sample oscillates spatially around 
a constant value. The oscillation vanishes in the non- 
interacting limit U — > 0. When g > (c^/i)" 1 , the general 
solution for the magnetic field becomes a superposition 
of the homogeneous and damping parts. We find that the 
damping part is energetically favored, thus, the Meissner 



effect occurs. Then, we consider the contribution from 
the helical edge state. As a result of Zeeman coupling, it 
is shown that the state exhibits weak diamagnetism that 
is independent of g. 

The argument we will present is focused on the quan- 
tized SH phase. Here, we mention that we have an upper 
limit for g(oc Ua 2 d). It has been pointed out that the 
system shows a phase transition from the quantized SH 
phase to the topological Mott insulating phase when U 
becomes larger. 5, 6,n ' The lattice constant a also should 
not be too large, since we consider the long-wavelength 
effective theory. 

This paper is organized as follows. In § 2, we introduce 
the layered KM model with on-site Hubbard U. In § 3, 
we integrate out Fermion and obtain a one-loop effective 
Lagrangian in the quantized SH phase. In § 4, we discuss 
the magnetic response. In § 5, we estimate the contri- 
bution from the helical edge state. In § 6, we comment 
on the relation to the superconductivity and our system. 
We use h = c = 1 unit and the Minkovskian metric 
g^ = diag{l, —1,-1), where //, v — 0, x, y. Summations 
run over repeated Greek indices. 

2. Layered KM Model with an Electron Corre- 
lation 

We consider electrons on the layered honeycomb lat- 
tice. We assume that interlayer coupling is negligibly 
small and each layer is described by the KM model. 6 ' 7 ' 
One of the essential ingredients of the KM model 1,2 ' is 
the intrinsic SO coupling Ago- We can say, not strictly 
but intuitively, that Aso gives an effective magnetic field 
depending on spin. It also gives an electronic excitation 
gap, 2 ' 



A = 3%/3A : 



so- 



(1) 
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Thus, as an analog of the quantum Hall effect, we see 
quantization of the SHC, 1 ' 2 ' 

e A 



■xy 



2nd\A\' 



(2) 



where d is the interlayer distance. The model can also 
have the Rashba extrinsic SO coupling A/j, which breaks 
the inversion symmetry and is induced by an electric 
field perpendicular to the honeycomb lattice plane. The 
term also breaks the conservation of electron spin S z . 
Hereafter, we consider the case of Ar = 0. 

We add the on-site Coulomb repulsion U > 0. The 
Hamiltonian per layer is 



H = 



1 E c \ c i + * Aso 



E 



'ij^i^z^j 



E< 



-U > n^riii, 



(3) 



where Cj (c|) is the annihilation (creation) operator of 
an electron with spin at the i-th site and t is the nearest 
neighbor hopping. The second term is the intrinsic SO 
term consisting of the next nearest neighbor hopping, 
and Vij = -js{di X d2) z = ±1, where di and d,2 are unit 
vectors along the two bonds where the electron moving 
from site j to i passes. 

Let us discuss how to deal with the electron correlation 
U. On-site Coulomb repulsion can be written by the on- 
site spin-spin interaction 



U , s U > t- ^2 

Um-frni = — (riit + n a) - Y^ c i SCt > ■ 



(4) 



The first term in the r.h.s merely gives the renor- 
malization for the chemical potential and can be ne- 
glected. We introduce the auxiliary field (p{, which is a 
three-component vector in the spin space, and use the 
Stratonovich- Hubbard transformation, 12 ' H — » Hsh = 
H + AH, where 



^ = ^(4^-^r, 



(5) 



H, 



sn 



t y^ cUi + «A S o 






E 



■'i3<-i°z^j 



«ij» 



E 



+ s 3 



^^r- 



(6) 



The spin-spin interaction is eliminated in appearance, 
but instead, we have coupling between (pi and the elec- 
tron spin, and a quadratic term of (pi. 

We consider the continuum limit and take into account 
the low-energy electronic excitations around K and K' 
points in the Brillouin Zone, 1,2 ' i.e., we omit the inter 
valley scattering. We introduce the electromagnetic U{\) 
gauge field A^ and SU(2) spin gauge field a M via the 
covariant derivative 



i-Dp = id^ - eA^ 



s 



(7) 



where do = Cp (the auxiliary field in the continuum limit) 
and a is an external field introduced artificially to esti- 



mate the spin current. We define a parameter 

4Ua 2 d 
9 = — z — ' 



(8) 



where a is the lattice constant, and the microscopic La 
grangian density is 13 ' 

£ = #+ {iD Q - iv(D x T z a x + D y a y ) + Ar z a z s z } * 



B 2 



1 



*, % WI - (9) 

where <3> = ^ Tas is the eight-component fermion field la- 
beled by the eigenvalues of the diagonal components of 
valley spin r, sublatticc spin a, and real spin s/2. The 
parameter v is the Fermi velocity when the system is in 
the metallic state, and eo and /zo denote the vacuum val- 
ues of the dielectric constant and magnetic permeability, 
respectively. Note that, except for the last term, the La- 
grangian (9) possesses the U(l) em x U(l) z local gauge 
symmetry. The SU(2) gauge symmetry is broken down 
to U(l) z , since the SO term contains s z . 

3. One-Loop Effective Lagrangian in the Quan- 
tized SH Phase 

The phase diagram for correlated electrons on pyro- 
clorc and honeycomb lattices with Aso has been dis- 
cussed using the slave- rotor model. 5 ' 6 ' The system with 
strong Aso an d small U is in the quantized SH phase 
when the spin conservation is preserved. As U increases, 
the band gap closes and the Mott gap opens instead, i.e., 
band-Mott transition occurs. We focus on the quantized 
SH phase. 

The derivation of the effective action in this phase is 
equivalent to that presented in ref. 14, although the phys- 
ical meaning of the spin gauge field is different. We inte- 
grate out \& from eq. (9) and obtain the one-loop effective 
Lagrangian for the gauge fields. We assume that the am- 
plitude of the gauge fields is small and use the Gaussian 
approximation. We consider fields with a length scale 
(temporal scale) of modulation that is sufficiently grad- 
ual compared with a (A), and use the long- wavelength 
(low-frequency) approximation. 

The result is 14 ' 
1 



£cfi — — 



25 



£ 



inch 



(10) 



Cb 



a xy 6 a u°P A v + 



eE 2 



2 2/x 2 2fx s 
+ (terms independent of a^ and A^), (11) 

where And stands for the induced part of the effective 
Lagrangian, e 012 = e 120 = e 201 = -e 021 = -e 210 = 
— e 102 = 1, and e s = — a z — Vag and b s = £)- e^did* 
with e 12 = — e 21 = 1 are the spin electric field and spin 
magnetic field, respectively. The first term in eq. (11) is 
the topological BF term, 7, 14 ~ 18 ) which plays an impor- 
tant role in our discussion. The coefficient is the quan- 
tized SHC given in eq. (2). Note that only a z couples 
to the electromagnetic gauge fields. This comes from the 
fact that the SU{2) symmetry is broken down to U(i) z 
symmetry by the SO coupling. Wc comment on the more 
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detailed properties of this term in the Appendix. The 
Maxwell term is renormalized as 14 ' 19 - ) 



(12) 
(13) 
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= £o + de, 
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6e e 


e 2 


67r|A|d' 
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e v 
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67r|A|d' 


e s = 


Se 
4e*' 


1 


1 


Ms 


4e 2 <5/i' 



By using the relations 



(14) 
(15) 

(16) 
(17) 
1/137 and e Mo = 1 



and also the parameters in Table. I, which are relevant 
for the honeycomb-layered insulator Na2lrC>3, 4 ) we ob- 

tain e o/e^W = °- 5 and Mo ■ ^£ld = 2 x 10_8 ' Le -> 
M ~ Mo- The elastic term for ag(= <p z ) is also induced. 
We can recognize that any potential terms (i.e., zeroth- 
order terms with respect to the derivative <9 M ) of A^ and 
also afj in And are absent because of the presence of 
U(l)om x U(l) z gauge symmetry in the fermionic part 
of the microscopic Lagrangian (9). Thus, the low-energy 
and long-wavelength physics of A^ and a z is described 
definitely by eq. (10). 



Table I. Parameters used for estimations. These are typical val- 
ues for Na2lr03, 4 ) which is a honeycomb- layered insulator with 
Ago and electron correlation. 



A 


U 


d 


a 


V 


0.5eV 


0.5eV 


10A 


1CL4 


3 x 10 4 m/s 



4. Magnetic Response 

We consider the static magnetic response. Here, we 
set a = const. The equations of motion for spin chemical 
potential a 2 , (originally, this is the auxiliary field <p z ) and 
magnetic field B obtained from eq. (10) are 



n 2 2 

e.V a n 



lag 



<yB, 



M .£-- 



CijVjB 



E 



eij\/ja . 



(18) 



(19) 



The r.h.s. of eqs. (18) and (19) are the results from the 
BF term [see eqs. (A-3) and (A-4) in the Appendix], re- 
spectively. 

We consider a sample in i > 0, and apply the ho- 
mogeneous magnetic field Bo parallel to z-axis, which is 
perpendicular to the layers. Around the boundary x = 0, 
we have the helical edge mode. 1,2 ) The contribution from 
the edge mode is discussed in the next section. It will be 
shown that the edge mode gives a small correction via 
Zeeman coupling. 



Obviously, the fields depend on x only, and the general 
solution is 



Go 0*0 
B{x) 



9°, 



■ni 



V^ y M ^ jk 



IKOX _i_ 



l /<y" j e -ifc x ) ( 20 ) 



(21) 
(22) 



(23) 

where a, and f3± are arbitrary constants. The energy 
functional of the fields is 

f d 3 x \J_ B 2 + ^L ( y a z ) 2 + 1 2I (24) 



Aft 


= 


c.r s 


C V a ' 


c 


= 


<,y-> 


s 


= 


<lw > o, 



2n 2C 



2.9 



which will be used to determine the constants. Note that 
the BF term is absent, since the term does not consume 



energy 



13) 



4-1 Oscillation effect 

For s = o-^ 2 fig < 1, Aft is real. The Dirichlet-type 
boundary condition for B{x) at x = is 

B(0) Bo 

M Mo' 

We also assume that the fields are real and that the 
energy functional (24) per period 2ir/k along the x- 
dircction takes a minimum value. Then, we obtain 



(25) 



B(x) = B { 



/' 



1 X 

+ cos k x } , (26) 



Mo U + A 1 + X 



alix) = B^^ 



1 X 



X 



Mo [1+X <7 s xy [il + X 
2s 2 (l + 2s) 



cosAftx X27) 



(28) 



1 + s 2 ' 

The first term in the r.h.s. of eq. (26) is homogeneous, 
and the second one shows oscillation of the magnetic 
field. It is reasonable that the amplitude of oscillation 
vanishes in the non- interacting limit U — > [see eqs. (8), 
(23), and (28)]. 

Using the parameters shown in Table I, we obtain s ~ 
7.0 x 10~ 6 , and it is hard to observe the oscillation. If we 
have the lattice with a = 350 nm, we obtain s ~ 0.99. In 
this case, the amplitude of oscillation is about 75 % of 
-Bom/Mo- The wavelength of the oscillation is 



2tt 



^osri. — ^ 



C V 1 - s 



3fim. 



(29) 



This result is consistent with a long-wavelength approx- 
imation since A osc i. ^ a. 
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4-2 Meissner effect 

For s > 1, feg becomes pure imaginary. We introduce 
a real value, 



«o = — ika = C 



1 



(30) 



We impose solution (20) to be real and finite, and use 
boundary condition (25). The solution that gives the 
minimum of energy functional (24) is 



B(x) -- 


= S — e- K0X , 
Mo 


(31) 


a z (x) = 


B e -K0X 


(32) 



^SyMO 

It is obvious that the energy of these solutions converges 
because of the exponential damping. The homogeneous 
part should not appear, since energy functional (24) di- 
verges. The above solutions remind us of the Meissner 
effect with the penetration depth 



A, 



2tt _ 2tt 

K C 



s-1 



(33) 



It has been pointed out that eqs. (18) and (19) lead 
to the London equation for s 3> 1." In this limit, we 
can neglect the 1/g term in the l.h.s. of eq. (18) and 
obtain the London equation by taking the rotation of 
both sides of eq. (19) and using eq. (18). Amazingly, we 
obtain the Meissner effect without the London equation 
in this paper. Namely, the condition for the Meissner 
effect is weakened as s > 1 , instead of s ^s> 1 . 

We can sec from eq. (33), large s strengthens the 
screening. On the other hand, U should not be too large, 
since we are discussing the quantized SH (topological 
band insulating) phase. 5, 6 ' Thus, we can see, from eqs. 
(2), (8), and (23), that large a and small d are favorable 
for observing the Meissner effect. We also note that large 
A oc Aso shortens the penetration depth, since C ex A 1 ' 2 
[see eqs. (14), (22), and (33)]. 

The physical picture of this Meissner effect is consid- 
ered to be as follows: the spin-orbit coupling opens the 
topological gap and compensates the energy loss coming 
from the screening of the magnetic field. 

If we have a lattice with a = 380 nm, instead of a = 10 
A in Table I, we obtain s ~ 1.012 and the penetration 
depth of the magnetic field is estimated to be A pcn . ~ 
2.8 fi m. This result is consistent with long- wavelength 
approximation, since A pC n. 3> a. 

4-3 Quantized current in topologically non trivial insu- 
lators 
The existence of charged current is indicated, since the 
magnetic field is spatially dependent. It may sound curi- 
ous that current flows in an insulating system. We note 
that the current is a result of a combination of spin ac- 
cumulation (A-3) and the dual quantized SH effect (A-4) 
derived from the BF term [see Appendix] . The spin accu- 
mulation causes the spin chemical potential a§, and the 
quantized electric current flows perpendicular to — Vag, 
i.e., the spin electric field. Therefore, the origin of the cur- 
rent has a direct analog of the quantized Hall effect, 20 ' in 



which non dissipative quantized transport carried not by 
the excited state but by the ground state occurs. These 
are typical transport phenomena in band insulators with 
non trivial topology 20,21 ' We note that such a non dis- 
sipative charge transport does not need the spontaneous 
U(l) em symmetry breaking. 

In this section, we have considered the bulk state only 
and the helical edge state was not taken into account. We 
will discuss, in the next section, that, as a result of the 
Zeeman effect, the edge state shows weak diamagnetism 
independent of the correlation parameter g, and gives a 
correction to the results of eqs. (26), (27), (31), and (32). 

5. Contribution from the Helical Edge State 

In this section, we consider the edge in more detail. 
It is important to estimate the contribution from the 
helical edge mode, which is the hallmark of the quan- 
tum SH system. 1, 2 > 8 ~ 10 ) We assume that the edge modes 
are completely localized at the boundary and neglect the 
broadness. First, we neglect Zccman coupling, and later, 
we take it into account. 

5.1 Without Zeeman coupling 

In this subsection, we omit Zeeman coupling and show 
that the helical edge mode has a non-gauge-invariant 
contribution to the magnetic response. This contribution 
is cancelled out by the anomalous boundary response 
as a result of the bulk BF term. Namely, only the bulk 
contribution taken into account in the previous discus- 
sion is relevant. This anomaly cancellation, renowned as 
the bulk-edge correspondence, is a natural consequence 
of the gauge invariance, and exactly the reason for the 
presence of the helical edge mode. 15-17, 22 ~ 26 ) 

We consider the layered system in i > 0, and each 
layer is normal to the z-axis. The helical edge mode 
can be modeled by a pair of quasi one-dimensional (ID) 
massless Dirac fermions with opposite spin and velocity. 
For simplicity, we assume that edge fermions are local- 
ized completely at the boundary. We also assume that 
the many-body interaction between edge fermions is ir- 
relevant, since we are discussing the integral quantized 
phase. We introduce gauge couplings for the fermions in 
a covariant manner, and integrate out the fermions. This 
calculation is a direct analog of that for the chiral edge 
mode in the quantized Hall effect, and can be modeled by 
a single ID massless Dirac fermion. 24 ' Then, we obtain 
the response to the electromagnetic gauge field 



■(edge)z 



<xf-f X 



(34) 



{ga' g/3' 
gap ~ (SW + e a a') g2 (gfl'fl ~ £ P' 73 ) f ^ \ 

where a,/3 = 0,y denotes the space-time coordinate at 
the quasi- ID edge lying in the x = plane, g a/3 = 
dtag(v- 2 ,-l), 3 2 = g a pd a dP = (d /v) 2 - d 2 yl and 
e Q/ g = — epa. This current is anomalous since it is not in- 
variant under the gauge transformation A a — > A a + <9 a £, 
where £ is a regular function. 

In the bulk, we have the BF term. The term in the 



J. Phys. Soc. Jpn. 



Full Paper 



Author Name 



symmetrized form can be written as 



>~BF — —pr 



0(x)e^(a^d p A u+ A p d p a z u ). (35) 



We see that the response to the external field A p has a 
boundary term [see the last term]: 

3C B F 



j{BF)z 



da z ^ 



(36) 



{x)a s xv e lxpv d"A"-8{x)- 






Note that this is spin density (/i = 0) or spin current 
density (/i = i). and the boundary term breaks the gauge 
invariancc. 

We consider the response to the static magnetic field. 
Around the edge x = 0, we may write Aq = A x = 
and A y ~ B(Q)x. Thus, A a = d a {B(0)xy} and from eq. 
(34), we obtain 



.(cdgc)z 
JO 

•(edgc)z 



S(x)^A y , 



(37) 
(38) 



and from cq. (36), 



jf F) * = e{x)al y B-5{x) a -fA yi (39) 



,-(BF) 2 



0. 



(40) 



Therefore, the non-gauge-invariant boundary terms can- 
cell each other out. Total spin density, which should be 
inserted in the r.h.s. of eq. (18), is 



.(BF)z .(edge)* 
J0 "T Jo 



(x)a s xy B. (41) 

Thus, we conclude that only the bulk contribution is rel- 
evant to the magnetic response when we neglect Zeeman 
coupling. 

5.2 With Zeeman coupling 

We can see that the helical edge mode contributes to 
the magnetic response as a result of Zeeman coupling. 21 > 
We assume that the Zeeman splitting energy is much less 
than the bulk band gap A. As mentioned, the helical edge 
mode is described by a pair of quasi-lD gaplcss Dirac 
fcrmions. When SHC is given by eq. (2), i.e., quantized 
as +1 with the unit of e/2-7r, the edge spectrum is 



E 



t,i 



Tvky, 



(42) 



where k y is the momentum measured from the Fermi 
points, and —v and +v (v > 0) are velocities for up-spin 
and down-spin fcrmions, respectively. In the case of SHC 
quantized as —1, the signs of velocities are opposite to 
each other. 

As usual metals, Fermi points for up- and down-spin 
fcrmions are split by the Zeeman effect, and the number 
density for each spin is generated as 



An t4 






(43) 



where [is is the Bohr magneton, and L (~^> a) is the total 
length of the sample edge. Induced spin density is 



A Ps 



-(An t -An 1 ) 



PbB{0) 
lixvd 



S(x). 



(44) 



Because of the characteristic feature of the spectrum 
(42), we have induced current density. For each spin, 

Aj 1 ^ = ±evAn t ' i , (45) 

and in total, we have 

Aj A/ + Aj l 

(46) 



ej fB(0)6(x) 
ira 



Therefore, the equations of motion (18) and (19) are 
modified as 






9 

1 dB 
[i dx 



y B B(Q) 

2nv 



6(x), (47) 



da^ efiBB(0) 



S(x). (48) 



xy dx nd 

Instead of eq. (25), the Dirichlet boundary condition for 
B at x = is 

-B(O) Bo ep B 

Mo 



-B(0), 



(49) 



/i no nd 

i.e., the value B(0) is shifted. Thus, the corrected results 
are given by replacing the coefficient B in eqs. (26), (27), 
(31), and (32) as 

Bo 



Bo 



1 



(50) 



(eii/iB/nd) 
We note that efXfXB/^d ~ 5.6 x 10~ 6 for /i = /io and 

d = 10A. 

The result indicates that the helical edge mode shows 
the weak diamagnctism independent of the parameter g. 

6. The Electric Conductivity: Comparision with 
Superconductivity 

Finally, we examine the electric conductivity in our 
system and compare it with the superconductivity. We 
integrate out do from the effective Lagrangian (10). Note 
that this integration is justified since do was originally 
introduced as the auxiliary field of the Stratonovich- 
Hubbard transformation. We do not integrate out a, 
since this field was introduced merely to estimate the 
spin current. This integration is straightforward, since 
the Lagrangian (10) is quadratic with respect to So- We 
obtain 

u s2 V7 2 

£eff - °x V 2^ e « ffl i L ~ ~Y A £sV 2 + -1 A 

ij 

eE 2 B 2 

+ ~ ~ 2^ 

+ (terms independent of af and A p ), (51) 

where Aj = Y^A^ij ~ dify /^ 2 ]Aj is the transverse (i.e., 
gauge invariant) component of A;. We set a, to be con- 
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stant after the integration. The quantized SH current is 
obtained from the first term, as expected. The electric 
current is 28 -* 



<lv 2 



e,V 5 



(52) 



We may take the gauge A — V • A = then we have 
E = — A and A T = A. Thus the electric conductivity is 



a xx (uj,q) 



S 2 2 

'xyl 



1 



e s q 2 -g 1 



(53) 



which vanishes in the DC limit taking q — > first 
and lo — > later when 1/g =^ 0. 12 * 1 Namely, the sys- 
tem is insulating. Then, we conclude that the system 
shows the Meissner effect without infinite DC conduc- 
tivity when g > (o x \ v fx)~ , in contrast to the super- 
conductivity. A similar result has been obtained by the 
Maxwell-Chern-Simons (MCS) theory, which is the low- 
energy and long- wavelength effective theory for the time- 
reversal- violating topological band insulator, i.e., the 
quantized Hall system. 29 ~ 32 ) 

Eqs. (52) and (53) indicate that the system becomes 
superconducting at the point 1/g = [see also refs. 14, 
33]. Actually, we can see the infinite DC conductivity 
from the Kramers-Kronig relation as Ke[a xx (uj, 0)] oc 
5(uj). Unfortunately, this point is difficult to realize in 
the quantized SH phase, as mentioned in Introduction. 

7. Summary 

In this study, we investigated the magnetic response in 
the quantized SH phase of a layered-honeycomb lattice 
system with the intrinsic spin-orbit coupling Aso and 
on-site Hubbard U. When g = AUa 2 d/3 < (cr^)" 1 , 
where a and d are the lattice constant and interlaycr dis- 
tance, respectively, the magnetic field inside the sample 
oscillates spatially around a constant value. The oscilla- 
tion vanishes in the non-interacting limit U —> 0. When 
g > (^x^/-*) -1 ' tne Meissner effect occurs. It may be pos- 
sible to see the oscillation or Meissner effect in a super- 
lattice system with an appropriately large a. As a result 
of Zeeman coupling, the helical edge state 1, 2 ) shows the 
weak diamagnetism that is independent of g. 
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Appendix: Physical implications of BF term 

In this Appendix, we summarized the properties of the 
BF term 7 - 14 - 18 ' in eq. (10) 

£bf = <r 8 xy e^ aii d p A v , (A-l) 

where a* is the quantized SHC given by eq. (2). 



A.l Quantized SH effect 

The spin current density obtained from the term is 

<9£bf 



J; 



daf 



<y Y, ei i E i- 



(A-2) 



3=x,y 

This shows the quantized SH effect. 

A. 2 Spin accumulation 
The spin density is 

<9£ B f 



daf, 



<4yB, 



(A3) 



where B is the magnetic field perpendicular to the honey- 
comb lattice layers. It resembles the Zeeman effect, but 
an essential difference is that the coefficient is not the 
Bohr magneton but the quantized SHC. 

A. 3 Dual quantized SH effect 
The electric current density is 

<9£bf 



dAi 



a xy X! e y' V J a 0' 



3=x,y 



(A4) 



This shows that the current flows perpendicular to the 
gradient of spin chemical potential ag, namely, the spin 
electric field. This may be called the dual quantized SH 
effect. 
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